Unconventional Superconducting Symmetry in a Checkerboard Antiferromagnet 
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We use a renormalized mean field theory to study the Gutzwiller projected BCS states of the 
extended Hubbard model in the large U limit, or the t-t'-J-J' model on a two-dimensional checker- 
board lattice. At small t' /t, the frustration due to the diagonal terms of t and J' does not alter the 
d x 2_ y 2-w&ve pairing symmetry, and the negative (positive) t' /t enhances (suppresses) the pairing 
order parameter. At large t /t, the ground state has a s — s wave symmetry. At the intermediate 
t' It, the ground state is d + id or d + is- wave with time reversal svmmetrv broken. 

PACS numbers: 74.20.Rp, 74.20.-z, 74.25.Dw 



I. INTRODUCTION 

Geometrically frustrated systems with strong correla- 
tion have attracted much attention due to the highly non- 
trivial interplay between frustration and correlation 1 ^. 
In such systems the pairwise interaction does not coin- 
cide with the geometry of the lattice, which may lead 
to exotic ground states. In particular, frustration in 
quantum magnets may cause certain types of magneti- 
cally disordered quantum phases, including the resonat- 
ing valence bond (RVB) spin liquid state 3 - and the valence 
bond crystal state 4 -. The quantum spin liquid state could 
become unconventional superconducting state when the 
charge carriers are introduced. There have been ex- 
perimental evidences for the unconventional supercon- 
ductors in these systems. Examples are the triangu- 
lar layer cobaltates compound Na^CoO^, layered or- 
ganic conductor k-(ET) 2 Cu2(CN)3^, the Kagome com- 
pound SrCrsGa40ig- , and 3-dimensional (3D) beta type 
transition-metal pyrochlore material KOs20e^. 

To describe the interplay between frustration and cor- 
relation, we consider a t-t'-J-J' model on a 2D checker- 
board lattice, and analyze the possible superconducting 
pairing symmetry of the model. The checkerboard lat- 
tice is a frustrated one, and may be considered as a 2D 
projection of a 3D corner-sharing lattice of pyrochlore. 
A schematic checkerboard lattice is illustrated in Fig. 1. 
The Hamiltonian reads, 

H = ~ £ hjPDiAaCja + h.C.)P D + ^ Jij Si ' Sj 

~ V^c\a C ™ (!) 
ia 

where c| CT is an electron creation operator with spin a at 
site i, Si is a spin operator for electron, /i is the chemi- 
cal potential, and (ij) denotes a neighboring pair on the 
lattice. Pd is a Gutzwiller projection operator to impose 
no double electron occupation at any site on the lattice. 
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FIG. 1: (Color online). Schematic structure of 2D checker- 
board lattice. Solid (open) circles represent sublattice A (B). 
The hopping integral and spin-spin coupling are t and J for 
the nearest neighbor links (solid lines), and t' and J' for the 
diagonal or next nearest neighbor links (dashed lines), respec- 
tively. 



Uj and Jij stand for the hopping integrals and antifer- 
romagnetic exchange couplings respectively, and iy = t, 
Jij = J for the nearest neighbor (n.n) links and = t' , 
J L j = J' for the diagonal or the next n.n. links as shown 
in Fig. 1. For convenience, we use x,y to represent the 
n.n. links while a, b to describe the two diagonal links. 
The Hamiltonian may be viewed as a strong coupling 
limit of a Hubbard model on the lattice with n.n. and 
next n.n. hopping integrals t and t' respectively and a 
on-site Coulomb repulsive interaction U . Hereafter we 
use t as an energy unit and set Jjt = 1/3. We choose 
J'/ J = (t' /t) 2 , consistent with the superexchange rela- 
tion of J = 4t 2 /U in the large U limit of the Hubbard 
model. 

The model has certain limiting cases. At t' jt — > 0, the 
model is reduced to the t-J model on a square lattice. At 
t'/t — ► oo, the model becomes a collection of independent 
ID t'-J' chains. At t'/t = 1, it is an isotropic checker- 
board lattice with a highly geometrically frustrated struc- 
ture. 



2 



Previous theoretical investigations mainly focused on 
the half-filled case without charge carriers. A variety of 
techniques have been employed to study the quantum 
antifcrromagetism on such a lattic o 10 ' 1:L i 12 ' 13 i 14 . Various 
quantum paramagnetic ground states may appear, which 
include some translational symmetry breaking states and 
a quantum spin liquid state. The introduction of dop- 
ing with mobile charge carriers in a frustrated quan- 
tum antiferromagnet may result in the appearance of 
unconventional superconductivit y 15 ' 16 i 17 ' 18 i 19 . Generally 
speaking, geometric frustration may play a key role in 
the mechanism of unconventional superconductivity. Re- 
cently exact diagonalization approaches have been em- 
ployed to study the superconducting fluctuations in this 
system with t'—t and J'= J 20 i 21 i 22 . They found evidence 
of enhancement in pairing amplitude at arbitrarily small 
J/t for a specific sign of the hopping amplitude. 

In this paper, we apply the plain vanilla version of the 
RVB theor y 23 i 24 to study the ground state of the t-t'-J- 
J' model on a 2D checkerboard lattice. The competition 
among various superconducting states will be examined. 
Since our primary interest is on the possible pairing sym- 
metry of the superconducting state for the doped system, 
we will not consider the possible long-range magnetic or- 
dering in the present paper. Our main results can be 
summarized below. The d x 2_ y 2-w&ve pairing found for 
the t-J model extends to a large region in parameter 
space of t' /t and doping concentration 8, and the nega- 
tive t' /t enhances the pairing while the positive t' /t sup- 
presses the pairing. At small doping and for \t' jt\ ~ 1, 
there is a small region where the pairing symmetry is 
d + id or d + is. At \t'/t\ > 1, there are regions where 
the pairing symmetry belongs to a s — s wave. 

The rest of the paper is organized as follows. In Sec. 
II, we apply the renormalized mean field theory to study 
the RVB state. In Sec. Ill, we present our numerical 
results on the possible superconducting ground states. 
In particular, four distinct superconducting phases show 
up in the phase diagram as a function of t' /t and the 
doping 8. Sec. IV is a summary. The diagonalization of 
the mean-field Hamiltonian and the explicit form of the 
self-consistent equations are presented in the Appendix. 



(1). In the renormalized mean field theory, we adopt the 
Gutzwillcr approximation to replace the projection by 
a set of renormalized factors, which are determined by 
statistical counting o 25 i 26 . We have 

( c L c j<?) = 9t{4 a c ja ) , (Si ■ Sj) = g s (Si ■ Sj) (3) 

where gt and g s are the Gutzwiller renormalized factors 
and are given by^ 3 - g t — 25/(1 + 5) and g s = 4/(1 + 8) 2 
where 8 denotes the doping density. Therefore, the vari- 
ational calculations of H of Eq. (1) in the projected BCS 
state is reduced to the variational calculations of an ef- 
fective Hamiltonian H e ff given below in the unprojected 
BCS states \^>bcs) for H eff . 



H, 



(Kl) 



(4) 



To proceed further, we introduce particle-particle and 
particle-hole mean fields, (r = ±x, ±y, ±a, ±6, see Fig. 
1) 



A T 



\ c i1; c i+Tl c iI c i+t]/0 

X^ c t Ci +T,<x)o- (5) 



Here we focus on the translational invariant state with 
the spin singlet and even parity superconducting pairing 
symmetry, where A M+T = A i+T:i = A r , and Xi,i+r = 
Xt- The superconducting order parameter is a 2 x 2 ma- 
trix representing the two sublattices, 



A sc (/c) 



(ct .cK 



. (ct c f - ) 

klA -klA> v k]A -klB' 

(ct c f - ) (ct C f - ) 

x k]B -k{A' x k]B -k[B' 



(0) 



Within the Gutzwiller approximation, it is related to A T 
by 



v , ?1 f ^a^osk a /\ x cosk x + AyCOSky i ■ _ 

sc{ > ~ 9t \ A x cosk x + A y cosk y A b cosk b l [ > 



II. FORMALISM 

We use a Gutzwiller projected BCS state 3 - as a trial 
wavefunction to study the ground state and the corre- 
sponding pairing symmetry of the Hamiltonian (1). The 
trial wavefucntion is of the form, 



I* 



GS) 



BCS) 



(2) 



where rii tCr — c] a Ci i<T , and the projection operator J^[(l — 

i 

ritfriii) removes the doubly occupied electron states on 
every lattice site i. We use a renormalized mean field 
theory 23 i 24 to calculate the energy of the Hamiltonian 



In the limit of 8 = 0, consistent with the fact that the 
ground state is a Mott insulator at half filling. At small 
8, we have Asc(k) ~ 8. 

In terms of these mean fields, the effective Hamiltonian 
can be expressed as 



H 



\l I? = Y ~^9sJij [£,i]c\a c J<? + AijCurCjt + h.C.) 



-gtUj {c% a C i° + h-c) - [1 Yj n i° + const - ( 8 ) 



There are eight independent complex mean-field 
parameters: £ x , £ y , A x , A y on the n.n. links and 
£ a , A a , A;, on the next n.n. links. Here we assume 
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FIG. 2: Phase diagram of t-t'-J-J' model on a checkerboard 
lattice in parameter space of hole density 8 and t'/t obtained 
by the renormalized mean field theory. The phases are indi- 
cated by their superconducting pairing symmetry defined in 
Table I, and J'/J = (t'/t) 2 . 

all the particle-hole mean fields £ to be real. We denote 
&x,y, 9 a ,b, and 9 Xt a as the relative phases of (A X: A y ), 
(A a ,A b ), and (A x ,A a ), respectively. The energy per 
site can be expressed in terms of the mean fields and is 
given by 

E gs = -2g t t(Z x + S v ) - gtt'ita + &) 

-l 9 sJ(e x +^ + \A x \ 2 + \A y \ 2 ) 

-4 9sJ ' ( ^ 2 + ^ + |Aa|2 + |Ab|2) - (9) 

The mean field parameters £ and A and the chemical po- 
tential fi can be determined by solving the self-consistent 
Eq. (5) together with an equation for the hole density. 
The mean field state at zero temperature can be obtained 
by the diagonalization of Hmf- We then determine the 
lowest energy state for each set of parameters t'/t and 
5. The detailed formalism of the diagonalization of Hmf 
and the self consistent equations can be found in Ap- 
pendix. These equations can be solved numerically, and 
the results are given in the next section. 

III. NUMERICAL RESULTS 

In this section, we present the numerical results of the 
self-consistent renormalized mean field theory for Hamil- 
tonian Eq. (1) on the checkerboard lattice. We will first 
discuss the phase diagram, then provide detailed analyses 
of the mean field parameters as functions of 8 for several 
typical values of t'/t. The phase diagram is shown in Fig. 
2. The ground state at 8 = is a Mott insulator. At finite 
doping, the ground state is a superconducting state, with 
four different types of pairing symmetry, as illustrated in 



TABLE I: The pairing symmetries (shown in Fig. 2 and 
used in the text) of the ground states of t-t'-J-J' model on a 
checkerboard lattice and their corresponding mean fields A T 
of Eq. (5). Q T , T i are the relative phase between A T i and A T 



Pairing symmetry 


mean field parameters 


d 


Az = -A y , A a = A 6 = 


d + id 


A x = -A y , A a = -A b , 9 x , a « tt/2 


s — s 


A x = Ay, A a = Ab, 9 x , a — 7T 


d + is 


A* = -Ay, A a = A 6 , e x , a « tt/2 



Table 1. Here we classify the pairing symmetry in terms 
of the relative phases between A y and A^, between A& 
and A a , and between A a and A x . Such a classification is 
consistent with the four-fold rotational symmetry in the 
Bravis lattice of the checkerboard structure. 

At the limit t'=J'=0, the model is reduced to the t-J 
model, and the ground state has a d x 2_ y 2 or ci-wave sym- 
metr y 23 i 27 at finite doping. This pairing state is robust 
against the next n.n. terms. As we can see from Fig. 2, 
the d-wave pairing symmetry of the ground state extends 
to a large region of both positive and negative values of 
t'/t. In such state, A a = A/, = 0, but there is an ad- 
ditional self-energy term arising from the next n.n. spin 
coupling. There are nodal quasiparticles, whose position 
are determined by the crossing of the lines cos k x = cos k y 
and the Fermi surface, similar to those obtained in the 
t-J model. Note that the d-wave pairing has been previ- 
ously found to be stable against the weak frustrations as 
studied by various author a 15 i 16 ' 19 on anisotropic triangu- 
lar lattices. At large \t' /t\, the ground state has a s — s 
wave pairing symmetry with A^ = A y and A a = A5. 
In that state, the relative phase between A a and A^ is 
7r. Between the above two regions, there is a small pa- 
rameter region around \t' /t\ = 1, where the ground state 
has a d + id pairing symmetry at small 8, and a d + is 
phase at larger S for negative t'/t. In both d + id and 
d + is states, the relative phases between A a and A x 
are close to 7r/2, and are weakly dependent of S, and 
the time reversal symmetry is spontaneously broken. We 
note that Hamiltonian (1) is asymmetric with respect to 
positive and negative values of t'/t, which is reflected in 
our phase diagram. Similar to the case in the t-J model, 
the particle-particle mean field amplitude A disappears 
at large 8, and the ground state becomes a normal metal. 
The details of the phase boundaries between the super- 
conducting and normal metallic states will be elaborated 
below. We also note that in the limit \t'/t\ — > 00, the 
ground state becomes ID-like. 

In Fig. 3, we display the amplitudes of the mean 
field parameters as functions of 8 at the symmetric point 
t'/t =1. At very low hole density, the ground state 
has d + id-wave pairing symmetry, and A a = i|A a |, but 
|A a | <C A x . As 8 increases, the ground state becomes a 
d-wave, and A a = A& vanishes. In comparison with the 
mean field amplitudes found in the t-J model, £, x is sim- 
ilar and insensitive to 8, but A^ decreases more rapidly 



FIG. 3: (Color online). Amplitudes of mean fields £ and A 
as functions of hole density 8 for parameters t = t' = 1, J = 
J' = 1/3. The ground state has d + id pairing (see Table 1) 
at < 5 < 0.02, and d x - y 2-w&ve pairing at 0.02 < 5 < 0.15, 
and is a normal state at 5 > 0.15. 
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FIG. 4: (Color online). Amplitudes of the mean fields £ and 
A as functions of 8 for parameters t = 1, t'/t = 1.1, J = 
1/3, J'/ J = 1.21. The ground state is a s — s wave at 0.02 < 
8 < 0.12, and a normal state otherwise. 



in the checkerboard model as 8 increases. The latter may 
be understood as the consequence of the non-zero value 
of £ Q in the present case, which increases rapidly as 8 in- 
creases. Note that the s — s wave state has a very close 
energy, although it is slightly higher than either d+id- or 
d- wave states. Recently, Poilblano^i performed a finite 
cluster exact diagonalization study of the t-J model on a 
checkerboard lattice. Some exotic states for positive t are 
found to have d x 2- y 2-, s — s- symmetries, which appears 
to be consistent with our results. 

The mean field amplitudes as functions of 8 for the 
model at t'/t = 1.1 are depicted in Fig. 4. The ground 
state is a s — s wave at 0.02 < 8 < 0.12. The pairing 



FIG. 5: (Color online). Amplitudes of mean fields £ and A as 
functions 8 in the case of t = 1, J = 1/3, t'/t = —0.5 J' /J = 
0.25. A a = A b = 0, and the ground state has a d-wave pairing 
symmetry(see Table 1). 



amplitudes disappear around S = 0.12, indicating the 
ground state to be a normal metallic state at 5 > 0.12. As 
we can see, the correlations along the diagonal directions 
(a and b) become more important than those along x 
and y directions at t'/t > 1. Note that the results near 
the half filled need to be cautious. At the half filling, 
the mean field ground state has only non-zero values of 
£ a and £/,, indicating that the state is a collection of the 
decoupled chains along the directions of a and b. This 
may attribute to the poor Gutzwiller approximation on 
ID systems^. 

In Fig. 5 and Fig. 6, we show typical c5 dependence 
of the mean field amplitudes for parameters t'/t < 0. In 
this case, £ Q and £t, have the opposite sign with £ x or 
£ y . At small value of \t' /t\, the ground state is a d-wave 
state, the amplitudes of the mean fields are illustrated 
in Fig. 5 for t'/t = -0.5. In that state, A a = A b = 
0, and £ a = are small and changes sign at S > 0.2. 
Interestingly, A^ decreases much slower as 8 increases, 
in comparison with that in the t-J model. This suggests 
that the superconducting state may extend to a much 
larger hole density. 

In Fig. 6, we plot the mean field amplitudes for t'/t = 
— 1. Away from half-filling, the amplitudes of mean fields 
change nonmonotonically and there exist three distinct 
pairing symmetries with respect to different doping lev- 
els. As hole density increases, the ground state evolves 
from the d + id-wave state (A tt = — A& = i|A Q |) to the 
d-wave state and to d + is- wave state (A a = A& = i|A a |), 
as we can see from the figure that the amplitude of A a 
decreases to zero around 8 = 0.02 and arises again at 
8 = 0.06. In contrast to the d + id-wave state, the am- 
plitude of A ((,) is comparable to A^ in the d + is-wave 
state. In comparison with the case of positive t'/t, we 
note that the suppression of A^ becomes much slower 
as 8 increases and the superconductivity appears more 
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FIG. 6: (Color online). Amplitudes of mean fields £ and A as 
functions 5 for the parameters t = 1, J = J' = 1/3, t' = — 1. 
The ground state is d + id-wave at < 8 < 0.02, d-wave at 
0.02 < 5 < 0.06, and d + is-wave at S > 0.06. 



favored for negative t'/t. This result is in agreement 
with previous studies for a t-J model on a triangular 
lattic e) 15 ' 16 . An intuitive physical understanding of such 
effect can be given as follows: for positive t'/t, the t' and t 
terms in kinetic energy match quite well, so that t' term 
may enhance the kinetic energy (make it lower), hence 
suppress the pairing amplitude; for negative t'/t, the t' 
term may introduce frustration in kinetic energy, hence 
enhances the pairing amplitude. The situation here is 
similar to that in cuprates, where the positive t'/t case 
corresponds to the electron-doped system while the neg- 
ative t'/t case corresponds to the hole-doped system. It 
is well known that the hole-doped system has a higher 
transition temperature while the electron-doped system 
has a lower one and a substantial region in doping with 



antiferromagentic long range order which we have not 
considered in the present paper for simplicity. 



IV. SUMMARY 

We have applied the renormalized mean field theory 
to study the Gutzwiller projected BCS ground state for 
the t-t'-J-J' model on a highly frustrated checkerboard 
lattice. As charge carriers are introduced, the half filled 
Mott insulator is evolved to resonating valence bond su- 
perconducting states, with four types of pairing symme- 
tries depending on the hole density and the parameter 
t'/t. We have found that the d x 2_ y 2 symmetry state pre- 
viously obtained for the t-J model is most stable in a 
large parameter region for small values of \t'/t\. In this 
region, the d- wave pairing order parameter is suppressed 
for positive t'/t and enhanced for negative t'/t. At large 
value of \tf /t\, the ground state is found to be so-called 
s — s wave. Around \t' /t\ — 1, we have found the time 
reversal symmetry broken states, d + id- and d + is- 
symmetries states are the most stable. The appearance 
of the exotic superconducting pairing symmetry in the 
t-t'-J-J' model on a checkerboard lattice may be viewed 
as the interplay between geometrically frustration and 
strong electron correlation. Our results may shed light 
on the understanding of the novel pairing symmetry fea- 
tures of highly frustrated systems. 

HXH thanks X. Wan for numerical assistance. YC 
would like to thank Z.D. Wang and Q.H. Wang for help- 
ful discussions. This research was supported by the RGC 
grants (HKU-3/05C, HKU7057/04P and HKU7012/06P) 
of Hong Kong SAR government, seed funding grant 
from the University of Hong Kong, and NSF in China 
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V. APPENDIX 

In the appendix, we present detailed calculations of 
the mean field theory. We first diagonalize the mean 
field Harniltonian (6) on the checkerboard lattice. Wc 
carry out the Fourier transformation of electron operator 
in the two sublattices A and £>, 



B22 = -g s J'AbCOsk b 
3 

B12 = B 2 \ = -g s J(A x cos k x + A y cos k y ) 
with k T — k ■ f. The matrix can be diagonalized by 
employing an unitary transformation as 



-yjt ~? 

'kT' 'kT' -fci' -kl 



a lf b lf a -W b -H 



where the k dependence of uf and vf are implied. The 
matrix element u and v satisfy the following conditions: 



Civ = ^-^=a J ; (T exp(ik-ri), i^A (10) 



k 

CllJ =J2-/^ b k* c Mik-n), ieB (11) 

k 

The summation over k — (k x ,k y ) runs in the reduced 
Brillouin zone. N is the number of sites in each sub- 
lattice. By using a spinor representation, the effective 
Harniltonian can be written as 



H eff = Y.^l M k^k + C0TlSt 



(12) 



where r?t = (at , bt ,a ,-,,6 ,-,) and the matrix M? 

'fc v kf kf -kV -k.[' k 

reads 



= / A(k) - n B(k) 
k \ Bt(fc) -A(k)+ t i 



(13) 



where both A{k) and B(k) are 2x2 matrices whose ele- 
ments are given by 

An = ~(2g t t' + -g s J'£ a )cosk a 
3 

^22 = -(2g t t' + -g s J'£ b )cosk b 

A\2 = A 2 i = -2g t t(cosk x + cosk v ) 
3 

--g s J{£, x cosk x + Zycosky) 
3 

Bn = -g s J'A a cosk a 



i i 
E U>2* + V 1*V2 = E u l v 2 * - u l v 2 = 



In terms of it's and v's the self-consistent equations for 
eight mean field and the hole density are given by, 



= n E +(!-/)«*] cos fc 



k,n 



k.n 



& = ]vE[/K| 2 + (l-/)K| 2 ]cOs(fc5), 



k.n 



A ° = ^E( 1 - 2 /)«* cosfc - 



k.n 



k.n 



1 ~ S = ^E/K| 2 + (l-/)b: 



n|2 
2l i 



k,n 



where the band index n runs over 1 , 2 and / is the Fermi 
distribution function f(E n ) which is step function at zero 
temperature. The diagonalization of Mr and the self- 
consistent equations are carried out numerically to obtain 
the phase diagram and the results reported in the text. 



